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Abstract 

We study mechanical systems subject to constraint functions that can be dependent 
at some points and independent at the rest. Such systems are modelled by means 
of generalized codistributions. We discuss how the constraint force can transmit an 
impulse to the motion at the points of dependence and derive an explicit formula to 
obtain the "post-impact" momentum in terms of the "pre-impact" momentum. 
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1 Introduction 



Mechanical systems subjected to nonholonomic constraints have received a lot of attention 
in recent years in the literature of Geometric Mechanics (see [|, H, ||, |1^, |14|, |16|, ^7\, ^ 



19 , 28, 31] and references therein). Indeed, the dynamics of nonholonomic mechanics 
have been described from several approaches: Hamiltonian, Lagrangian and even Poisson 
methods have been used. 

The constraints which are usually considered in the literature (both linear and nonlinear) 
satisfy a certain regularity condition. That is, they are given by a set of independent nonholo- 
nomic constraint functions, or, in a global description, by a distribution on the configuration 
manifold in the linear case or a submanifold of its tangent bundle in the case of nonlinear 
constraints. 

However, there is an increasing interest in engineering and robotics in the motion of special 
mechanical systems as, for example, dynamical devices that locomote with the enviroment 
via impacts, sudden changes of phase space, etc. In many cases, the jump of the system's 
velocity is produced by an impulse that enforces new constraints on the system. In some 
cases, these systems admit a nice mathematical modelling. 

In this paper, we are interested in the following situation. In a local description, given a set of 
constraints . . . , ^m}, one assumes that they become linearly dependent at some points. 
In a global picture, the constraints are given by a generalized codistribution with variable 
rank. One could think of simple examples that exhibit this kind of behaviour. For instance, 
imagine a rolling ball on a surface which is rough on some parts but smooth on the rest. 
On the rough parts, it will roll without slipping and, hence, nonholonomic linear constraints 
will be present. However, when the sphere reaches a smooth part, these constraints will 
disappear. 

The first (and, up to our knowledge, the unique) reference in a geometrical context for 
such kind of systems is 0. In that paper, the constraints are provided by a set of global 
1-forms on the configuration manifold and, using the Frobenius theorem, the authors gave 
a classification of them according to the existence of some special sets that can exert a 
big influence on the trajectories of the system. In particular, the existence of an integral 
manifold gives a sort of partial holonomicity with strong implications. The authors were 
mainly motivated by problems in motion planning. However, we are interested, at least 
in this first approach, in the geometrical and topological aspects of the problem. In other 
words, we are concerned with obtaining the dynamical laws that govern the motion of the 
system. 

In consequence, we consider nonholonomic constraints given by a generalized codistribution 
D, that is, a codistribution which does not necessarily have the same rank at all points in the 
configuration manifold. This approach leads us to the definition of the concepts of regular 
and singular points. It should be noticed that our definitions are slightly different from 
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those in . Indeed, the regular points are those where the codistribution has locally constant 
rank. In this sense, the generalized codistribution is a regular codistribution, as is commonly 
understood, on the connected components of the set of regular points. The singular points 
are those where the codistribution changes its rank. From a dynamical perspective, the 
situation on the regular points is already known: we can derive the equations of motion 
following d'Alembert principle and treat them making use of the well-developed theory for 
nonholonomic Lagrangian systems. 

However, on the singular points the matter is essentially different. The classical derivation 
of the equations of motion no longer works and we must solve the problem with other 
methods. Here, we have adopted a point of view strongly inspired by the theory of impulsive 
mechanics ||, ||, |10|, |Tl], |2l|, and we use Newton's second law in its integral 



form [Q. Analyzing the trajectories which cross the singular set, we have found that, in 
certain cases, the constraint force can transmit an impulse to the motion. It is precisely 
the sudden appearence of new constraints (that is, the change of rank of the codistribution) 
which induces this impulsive character. More precisely, given a motion q{t) crossing the 
singular set at time to, we define two vector subspaces of T*^^^^Q as follows: is the 

limit of all the 1-forms in the codistribution based on q{t), t < Iq, when t ^ . D^^^^^ 
is defined analogously. Our conclusion is that there exists a jump of momentum only if 
is not contained in D~(^^^-^•^ and the "pre- impact" momentum p(to)- does not satisfy the 
constraints imposed by D'^i^t.o)' ^'^ch a case, we propose that the jump is determined by 
Apito) G D'^(^t^^ and the condition that the "post-impact" momentum p(to)+ must satisfy the 
constraints imposed by D~^{t^y 

To find out about the relation between the theory developed here and the Hamiltonian theory 
of impact [23, Q applied to this problem, which seems to be a promising possibility, is the 
object of current research. 

The paper is organized as follows. In Section 2, we introduce the notion of generalized codis- 
tribution, which is just the geometrization of constraints with non-constant rank. In Section 
3 we review the theory of impulsive forces and impulsive constraints. The ideas exposed 
here will be helpful in understanding the developments of Section 4, which constitutes the 
main contribution of this paper, where we study the equations of motion for mechanical sys- 
tems subjected to generahzed constraints and we derive the equations describing the jump 
of momenta. Finally, in Section 5, some examples are discussed with detail. 



2 Generalized codistributions 

We introduce here the notion of a generalized codistribution. This notion will be helpful 
in subsequent sections to model geometrically the dynamical systems under consideration, 
that is, systems subjected to constraints which can "degenerate" at certain points. All the 



results in this section are adapted from the ones stated for generalized distributions in [29 
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By a generalized codistribution we mean a family of linear subspaces D = {Dg} of the 
cotangent spaces T*Q. Such a codistribution is called differentiable if Vg G DomD, there is 
a finite number of differentiable local 1-forms cui, defined on some open neighbourhood 
U of q such that Dq/ = span{c<Ji(g'), ...,ui{q')} for all q' G U. 

We define the rank of D at g as p{q) = dim Dq. Given go G Q, if -D is differentiable, it is 
clear that p(g) > p(go) in a neighbourhood of go- Therefore, p is a lower semicontinuous 
function. If p is a constant function, then D is a codistribution in the usual sense. 

For a generalized differentiable codistribution D, a point q E Q will be called regular if g 
is a local maximum of p, that is, p is constant on an open neighbourhood of g. Otherwise, 
g will be called a singular point of D. The set R of the regular points of D is obviously 
open. But, in addition, it is dense, since if go G 5* = Q \ -R, and f/ is a neighbourhood of go, 
U necessarily contains regular points of D {p\u must have a maximum because it is integer 
valued and bounded). Consequently, go G R. 

Note that in general R will not be connected, as the following simple example shows: 

Example 2.1 Let us consider Q = and the general differentiable codistribution D(^x,y) = 
span{0(a;)((ix — dy)}, where is defined by 

(f){x) = { 1 ~ 
^ ^ I e"^ X > 

The singular points are those of the y-axis, and the connected components of R are the 
half-planes a; > (where the rank is 1) and a; < (where the rank is 0). 

Remark 2.2 We note that the notion of singular point defined here is different from the 
one considered in In that paper, the authors treat the case of generalized constraints 
given by a globally defined set of 1-forms, tui, tu;. Then, they consider the /-form 

VL = oji /\ ... /\ oJi . 

The singular set consists of the points for which fi(g) = 0, that is, the points, g, such that 
{u;i(g), u;;(g)} are linearly dependent. Applying this notion to the former example, the 
set of singular points would be the half-plane {x < 0}. 

Given a generalized codistribution, D, we define its annihilator, as the generalized 
distribution given by 

D° : DomDcQ — > TQ 

q ^ Dl = {Dqr. 

Remark that if D is differentiable, D° is not differentiable, even continuous, in general 
(the corresponding rank function of D° will not be lower semicontinuous). In fact, D" is 
differentiable if and only if D is a regular codistribution. 
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We will call M an integral submanifold of D if T^M is annihilated by Dm at each point 
m G M. M will be an integral submanifold of maximal dimension if 

TmM° = Dm, \fm e M . 

In particular, this implies that the rank of D is constant along M. A leaf L of D is a 
connected integral submanifold of maximal dimension such that every connected integral 
manifold of maximal dimension of D which intersects L is an open submanifold of L. D will 
be a partially integrable codistribution if for every regular point q & R, there exists one 
leaf passing through q. D will be a completely integrable codistribution if there exists a 
leaf passing through q, for every q & Q. In the latter case, the set of leaves defines a general 
foliation of Q. Obviously, any completely integrable codistribution is partially integrable. 

M being an integral submanifold of D is exactly the same as being an integral submanifold 
of its annihilator D°, and so on. 

In Example 2.1, the leaves of D are the half-plane {x < 0} and the half-lines of slope 1 
in the half-plane {x > 0}. Given any singular point, there is no leaf passing through it. 
Consequently, D is not a completely integrable codistribution, but it is partially integrable. 



3 Impulsive forces 

In this section, we discuss classical mechanical systems with impulsive forces [|I|, |ll|, ^ 



24| , This field has traditionally been studied by a rich variety of methods (analytical, 
numerical and experimental), being a meeting place among physicists, mechanical engineers 
and mathematicians (for an excellent overview on the subject, see 0]). Recently, such 
systems have been brought into the context of Geometric Mechanics [|, ^ |10|, |15[. We will 
give here a brief review of the classical approach. These ideas will be useful in understanding 
the behaviour of the constraint forces acting on mechanical systems subject to generalized 
constraints. Both situations are not the same, but have many points in common, as we will 
see in the following. 

Consider a system of n particles in such that the particle r has mass m,.. Introducing 
coordinates (g^''~^, g^''"^, g^'') for the particle r, we denote by Q the configuration manifold 
i?^" and by Fr = (F^''~^, F^''"^, F^^) the resultant of all forces acting on the r^^ particle. 

The motion of the particle r in an interval [t, t'] is determined by the system of integral 
equations 

mr{q'{t')-q'it))= f F\r)dr , (1) 



where 3r — 2 < /c < 3r and k is an integer. The integrals of the right-hand side are the 
components of the impulse of the force Fr. Equation @j establishes the relation between the 
impulse and the momentum change, i.e. "impulse is equal to momentum change". Equation 
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(H) is a generalized writing of Newton's second law, stated in integral form in order to allow 
us to consider the case of velocities with finite jump discontinuities. This is precisely the 
case of impulsive forces, which generate a finite non-zero impulse at some time instants. 

If F is impulsive there exists an instant to such that 

lim ! FMdr = P y^O. (2) 

Equation (H) implies that the impulsive force has an infinite magnitude at the point Iq, but 
we are assuming that its impulse P is well defined and bounded. The expression P ■ 5(to) = 
lim F{t) can be mathematically seen as a Dirac delta function concentrated at to- 

Now, we will derive the equations for impulsive motion following the discussion in |2^. In 



the sequel, the velocity vector of the r*'' particle, {q^'" ^,q'^^'), will be denoted by q^. 

Then, the system of integral equations can be written as 



m.(g^(to + e) - q'ito - e)) = T"^' F,(r)rfr . 

Jtn — e 



"to+t 
'to-e 

If we multiply this expression by the virtual displacements at the point g(to)) we obtain 

{Pr{to + e) - Prito - e)) ■ Sq"^ = / Fr{T)dT ■ dq' . 

Jto-e 

For the entire system, one has 

E Prito + e) - Pr{to - e) - / F;{r)dr \-Sq^ = J2 F';{r)dr ■ 6q^ , (3) 

where F^. and F" are, respectively, the resultant of the given forces and of the constraint 
reaction forces acting on the r*'* particle at time r. 

Now, take a local chart (q^), 1 < A < 3?7, on a neighbourhood U of qito) and consider the 
identification TgQ = i?^", which maps each Vg G T^Q to (v^), such that Vg = ^^"j^ ; 

each q E U. Let us suppose that the constraints are given on U by the 1-forms tOi = fiiAdq^, 
1 < i < m. Then, we have that fJ,iA{(l(t)) = f^iAilito)) + 0(t — to) along the trajectory q(t). 
As the virtual displacements at the point q{t) satisfy by definition 

J2l^iA{q{mSqit))^ = 0, l<t<m, 

we conclude that ^ (fJ'iAilito)) {6q(t))"^ + 0(t — to)) = 0. Therefore, we have that 

6q'{t) = 6q'{to) + 0(e) , t G [to - e, to + e] , 
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that is, the virtual displacements at q{t) can be approximated by the virtual displacements 
at q{to). As a consequence, in the right-hand side of (|) we can write 

/•to+e rtQ+^ A-to+e 

/ F;'{T)dT-6q' = / F;.'{T)-5q'dT= / F;'(r) ■ (5g"(r)rfr + 0(e) . 

J tQ — e Jto—€ Jto — e 

The first term after the last equality is the virtual work done by the constraint forces along 
the trajectory, and this work is zero since we are considering ideal constraints. The second 
one goes to zero as e tends to zero. 

In the presence of given impulsive forces acting on m particles, say, at time to, we have 

hm / Fr'{T)dT = P,, y^O , 1 <r' <m . 

J to 

Then, taking the limit e — in (|^), we obtain the equation for impulsive motion |2l|, 

n 

Yl {Prito)+ - PM- - Pr} ■5q'=Q. (4) 

r=l 

An example in which equation (^) can be applied is when we strike with a cue a billiard 
ball which is initially at rest. In that case we are exerting an impulsive force that puts the 
ball into motion. But what happens when the ball collides with the edge of the billiard? 
What we see is that it bounces, i.e. it suffers again a discontinuous jump in its velocity. 
The constraint imposed by the wall of the billiard exerts an impulsive force on the ball. 
When the impulsive force is caused by constraints, such constraints are called impulsive 
constraints. There is a number of different situations in which they can appear. In the 
following, we examine them. 

In the presence of linear constraints of type \& = 0, where ^ = hk{q)f (a situation which cov- 
ers the case of unilateral holonomic constraints, such as the impact against a wall, and more 
general types of constraints such as instantaneous nonholonomic constraints), the constraint 
force, F = Fk dq^, is given by 

Fk = H-bk , 

where is a Lagrange multiplier. Then the constraint is impulsive if and only if 

lim / jj-bkdr = 7^ , 

J to 

for some k. The impulsive force may be caused by different circumstances: the function 
is discontinuous at to; the Lagrange multiplier /i is discontinuous at to or both. 

The presence of such constraints does not invalidate equation (^). It merely means that the 
virtual displacements 6q^ must satisfy certain additional conditions, which are just those 
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imposed by the constraints. So, in the abscence of impulsive external forces and in the 
presence of impulsive constraints, we would have 



Y.^PM-Sq'^ = 0, (5) 



where Ap,.(to) = Pr(to)+ -Pr(to)-- 

Remark 3.1 In general, equation (^) is not enough to determine the jump of the momen- 
tum. One usually needs additional physical hypothesis, related with elasticity, plasticity, etc. 
to obtain the post-impact momentum. In this respect, there are two classical approaches, the 
Newtonian approach and the Poisson approach 0, ^ . The Newtonian approach relates the 



normal component of the rebound velocity to the normal component of the incident velocity 
by means of an experimentally determined coefficient of restitution e, where < e < 1. 
Poisson approach divides the impact into compression and decompression phases and relates 
the impulse in the restitution phase to the impulse in the compression phase. 



Remark 3.2 It could happen that impulsive constraints and impulsive forces to be present 
at the same time. For example, in the collision between a rigid lamina and an immobile plane 
surface, we must take into account not only the normal component of the contact force, but 
also the friction force associated to the contact. It is not innocuous the way the friction is 
entered into the picture. In fact, the Newton and Poisson approaches have been revealed to 
be physically inconsistent in certain situations. On the one hand, Newton approach can show 
energy gains [0, . On the other hand, Poisson's rule is not satisfactory since non-frictional 



dissipation does not vanish for perfectly elastic impacts 0, ^ . This surprising consequence 
of the impact laws is only present when the velocity along the impact surface (slip) stops 
or reverses during collision, due precisely to the friction. Stronge ^ proposed a new 
energetically consistent hypothesis for rigid body collisions with slip and friction. It should 
be noticed that the three approaches are equivalent if slip does not stop during collision and 
in the perfectly inelastic case (e = 0). 

Recently, a new Newton-style model of partly elastic impacts has been proposed ||2^ which. 



interestingly, always dissipates energy, unlike the classical formulation of the Newtonian 



approach discussed in p6|. 



In the frictionless case, one can prove the following 



Theorem 3.3 (Carnot's theorem) (/^^ WAD The energy change due to impulsive con- 
straints is always a loss of energy. 
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4 Mechanical systems subjected to generalized con- 
straints 



In this section, we study the equations of motion for mechanical systems subjected to general- 
ized constraints. Let us consider a mechanical system with Lagrangian function L : TQ R, 
L{v) = ^g{v, v) — {U o tq){v), where g is a Riemannian metric on Q and f/ is a function on 
the configuration space Q (the potential). Suppose, in addition, that the system is subjected 
to a set of constraints given by a generalized differentiable codistribution D on Q, that is, 
we assume that tq{D) = Q. The motions of the system are forced to take place satisfying 
the constraints imposed by D. 

We know that the codistribution D induces a decomposition of Q into regular and singular 
points. We write 

Q = RUS. 

Let us fix Rc, a connected component of R. We can consider the restriction of the codis- 
tribution to Rc, Dc = D\R^ : Rc C Q — > T*Q. Obviously, we have that Dc is a regular 
codistribution, that is, it has constant rank. 

Then, let us denote by D° : Rc — > TQ the annihilator of Dc. Now, we can consider the 
dynamical problem with regular Lagrangian L, subjected to the regular codistribution D° 



and apply the well-developed theory for nonholonomic Lagrangian systems [Q, 0, [T^, 

Consequently, our problem is solved on each connected component of R. The situation 
changes radically if the motion reaches a singular point. The rank of the constraint codis- 
tribution can vary suddenly and the classical derivation of the equations of motion for non- 
holonomic Lagrangian systems is no longer valid. Let us explore the behaviour of the system 
when such a thing occurs. 

Consider a trajectory of the system, q{t), which reaches a singular point at time to^ i-e. 
g(to) £ S, such that q{to — e, to) C R and g(to)^o + e) C -R for sufficiently small e > 0. The 
motion along the trajectory q{t) is governed by the following equation, which is, as in the 
impulsive case, an integral writing of Newton's second law, to consider possible finite jump 
discontinuities in the velocities (or the momenta). That is, at each component 

PA{t')-pA{t) = [ FA{r)dT, (6) 



on any interval t < t' < oo, where F is the resultant of all the forces action on the trajectory 
q{t). In our case, the unique forces acting are the constraint reaction forces. 

The nature of the force can become impulsive because of the change of rank of the codis- 
tribution D. We summarize the situations that can be found in Table 0. On entering the 
singular set, the rank of the codistribution D at the singular point g(to) can be the same as at 
the preceding points (Case 1) or can be lower (Cases 2 and 3). In these two latter situations. 
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the constraints have collapsed at qito) and this induces a finite jump in the constraint force. 
As the magnitude of the force is not infinite, there is no abrupt change in the momenta. 
Consequently, in all cases, we find no momentum jumps on entering the singular set. 





q{to — e): preceding points 


q{to): singular point 


q{to + e): posterior points 


Case 1 


p — r 


Po = ro = r 


p > r 


Case 2 


p — r 


Po^ro<r 




Case 3 


p — r 


Po = ro < r 


P>ro 



Table 1: Possible cases. The rank of D is denoted by p 



On leaving the singular set, the rank of D at the posterior points can be the same as at q{to) 
(Case 2) or can be higher (Cases 1 and 3). In Case 2 nothing special occurs. In Cases 1 and 
3, the trajectory must satisfy, immediately after the point q{to), additional constraints 
which were not present before. It is in this sense that we affirm that the constraint force 
can become impulsive: if the motion which passes through the singular set and tries to enter 
the regular one again does not satisfy the new constraints, then it experiences a jump of 
its momentum, due to the presence of the constraint force. In this way, the new values 
of the momentum satisfy the constraints. But one has to be careful: the impulsive force 
will act just on leaving S, on the regular set. Consequently, we must take into account 
the virtual displacements associated to the posterior regular points. The underlying idea 
of the mathematical derivation of the momentum jumps in Section 4.1 is the following: 
to take an infinitesimal posterior point q{t) to q{to), to forget for a moment the presence 
of the constraints on the path q{tQ,t) and to derive the momentum jump at q{t) due to 
the appearence of the additional constraints. Afterwards, to make a limit process t — > 
cancelling out the interval {to,t) where we "forgot" the constraints. In any case, we will 
make the convention that the jump happens at q{to). 

We illustrate the above discussion in the following example. 

Example 4.1 Consider a particle in the plane subjected to the constraints imposed by the 
generalized codistribution in Example 2.1. The Lagrangian function is 

L : TR^ — > R 

{x,y,x,y) I — > -m{x'^+y^). 

On the half-plane Ri = {x < 0} the codistribution is zero and the motion is free. Conse- 
quently the trajectories are 

X = Xot + Xo, 

y = ilot + yo- 
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If the particle starts its motion with initial conditions Xq = —1, i/q = 1, Xq = 1, yo = 0, 
after a time 1, it reaches the singular set S = {x = 0}. If the motion crosses the y-axis, 
something abrupt occurs on entering the half-plane R2 = {x > 0}, where the codistribution 
is no longer zero and, indeed, imposes the additional constraint x = y (Case 1). We know 
that the integral manifolds of D on R2 are half-lines of slope 1, so the particle suffers a finite 
jump in the velocity on going through the singular part in order to adapt its motion to the 
prescribed direction. 

If, on the contrary, the particle starts on R2, say with initial conditions xq — 1, yo — 1, 
xq — —1, yo — —1, after a certain time, it reaches the set S. On crossing it, nothing special 
happens, because the particle finds less contraints to fulfill, indeed, there are no constraints 
(Case 2). Its motion on Ri is free, on a straight line of slope 1 and with constant velocity 
equal to the one at the singular point of crossing. 



4.1 Momentum jumps 

Now, we derive a formula, strongly inspired by the theory of impulsive motion, for the 
momentum jumps which can occur due to the changes of rank of the codistribution D in 
Cases 1 and 3. 

At q{to) we define the following vector subspaces of T*^^^-^Q 



L>g(j^) = {a e T*^^^^Q / 3a : {t^ - e, t^) T*Q, a{t) e D^^t) and lim. a{t) = a} 




Prom the definition of D,. ^ and D„,. ^ we have that 

where denotes the orthogonal complement with respect to the bilinear form induced by 
the metric g on the cotangent space T*^^^^Q, and the hmits (-D"*")" and (-D"*")"*" are defined 
as in the case of D~ and . In the following, we shall not make a notational distinction 
between the metric g and the induced bilinear form on T*^^^^Q. In each case, the precise 
meaning should be clear from the context. 

Since D is a differentiable codistribution then 

Dqit,) C and D^ito) ^ D+^,) ■ 

Along the interval [to, t], we have 

PA{t) - pA{to) ^ [ FA{r)dr. 

J to 
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Multiplying by the virtual displacements at the point q{t), we obtain 

(pAit) - PAito)) ■ = /* ^^(^) ■ Kit) ■ (7) 

J to 

Summing in A, we get 

n 

Y^^PAit) 

A=l 

Since we are dealing with ideal constraints, the virtual work vanishes, that is 

A=l "^*o 

If t is near to, then r is close to t, and g(r) is near q{t), so 6q remains both nearly constant 
and nearly equal to its value at time t throughout the time interval {to,t\, in the same way 
we exposed in Section 3. Therefore, 

E / FAir) dr ■ 5ql^,^ = / ^A(r)5gt(,) dr + 0(t - to) = 0(t - to) • 

A=l ^"''O / yl=l -^to 

Consequently, equation (||) becomes 

n 

Y.(PAit) - PAito)) ■ Sq^^t) = 0{t - to) . (9) 

A=l 

Taking limits we obtain 




which implies 

n 

J2(PA{to)+ - PAito)) lim = (10) 

A=l *^*o 

or, in other words, 

iPAito)+ - PAito)) dq"" e lira D^^t^ = D+^^^ . (11) 



PAito)) ■ Sql^t) = H FAiT)dT ■ 5q 
A=l 



A 

m ■ 
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there is no jump of momenta 




possibility of jump of momenta 



Table 2: The two cases 

Conclusion: Following the above discussion, we will deduce the existence of jump of mo- 
menta depending on the relation between -0"^^^-, and D'^{t^^y The possible cases are shown in 
Table |. 

In the second case in Table we have a jump of momenta if the "pre-impact" momentum 
p{to)- = pito) does not satisfy the constraints imposed by D~^(t^y that is 

PA(t,)-dq^ i (DJ,^/ . 
Our proposal for the equations which determine the jump is then 

The first equation has been derived above (cf. (|TT])) from the generalized writing of New- 
ton's second law (|T|). The second equation simply encodes the fact that the "post-impact" 
momentum must satisfy the new constraints imposed by D~^^^^y 

Remark 4.2 In Cases 1 and 3, the virtual displacements at g(to) are radically different from 
the ones at the regular posterior points, because of the change of rank. From a dynamics 
point of view, these are the "main" ones, since it is on the regular set where an additional 
constraint reaction force acts. As we have seen, the momentum jump happens on just 
leaving S", due to the presence of this additional constraint force on the regular set. Note 
that with the procedure we have just derived, we are taking into account precisely the virtual 
displacements at the regular posterior points, and not those of g(to)- If we took the virtual 
displacements at q'(to) and multiply by them in (|^), we would obtain non-consistent jump 
conditions. This is easy to see, for instance, in Example |4.1| . 



An explicit derivation of the momentum jumps for Cases 1 and 3 would be as follows. Let 
m be the maximum between p = r, the rank at the regular preceding points, and p = s, 
the rank at the regular posterior points. Then there exists a neighbourhood U of q{to) and 
1-forms uJi, ...,u!m such that 

Dg = span {cji(g), ...,uJm{q)} , E U . 

Let us suppose that lji,...,Us are linearly independent at the regular posterior points (if 
not, we reorder them). Obviously, at q{to), these s 1-forms are linearly dependent. In the 



13 



following, we will denote by uoi the 1-form evaluated at (t time immediately posterior 
to to) i-e. = uJi{q{t)), in order to simplify notation. 

Since the Lagrangian is of the form L = T — U, where T is the kinetic energy of the 
Riemannian metric g, that is, L = IgABQ^Q^ ~ U{q), then we have that 



A{q{m\t) = Y,^^Ag^^pB{t) = 0, J = . (12) 



A.B 



Using the metric g we have the following decomposition 

t;q = d,®d^, qeQ. 

The two complementary projectors associated to this decomposition are 

V, : t;q 

Q, : T;Q — . D, 

The projector Vg is given by 

Vq{aq) =ag- C^aq{Zi)uj, ag G T*Q 

where 



Zi = g^'^uj. 



JiB 



d 



and are the entries of the inverse matrix of C, the symmetric matrix with entries Cj 
^iAg^^^jB-i or C = ujg~^uj'^ with the obvious notations. 

By definition 

PA{to)+dq\(^to) = lim(pA(t)rfg^|g(t)) • 



From ([T2|), Pg(t)(pyi(t)rfg^|g(t))) = PA{t)dq^\g^t) and then 

PA{to)+dq\(^to) = ^hm iPA{to)+dq\(^to)) ^ (^Jt^))^ • (13) 

Combining (pUD and (|T^) , we obtain 

PA(^o) + C?g^|g(to) = ^/^^'^eW^ [PA(io)-rfg^l9{to)] • 

In coordinates, this can be expressed as 



- Inn ( 

*^*o \i,j,A,B 



PA{to)+ = PAito)- - lim I C^ujBg^'^uJiA Pcito)^ , A = l,...,n. (14) 

/ lit) 
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Equation (p^ can be written in matrix form as follows 

p{to)+ = (^Id - \imiuj''C-'Lug-\^t)^ p(to)- ■ (15) 

With the derived jump rule, we are able to prove the following version of Carnot's theorem 
for generalized constraints. 

Theorem 4.3 The kinetic energy will only decrease by the application of the jump rule (fJ^j. 
Proof: We have that 

g {p{to) + ,p{to)+) = g \ \ lim Pg(t) p(to)-,p(to)- - hm Qg(t) p(to) 



g I I lim p{to)-,pito) 



) 



9 (p(^o)-,p(io)-) - ^ I I ^lim Qg{t) I p(io)-,p(ioJ 



Since 6f |^|^hm Qq(t) j p{to)-,p{to)- j = 5- |^ J^^lim Qq{t) j p{to)-, J^hm Qg(t) j Pito)-j > 0, we 

can conclude that ^ ^ 

2^(p(^o)+,p(^o)+) < -gip{to)-,p{tQ)-) . 

\qed\ 

In fact, the jump rule (0) has the following alternative interpretation. Let p G D^{tQ) ^^"^ 



observe that 

9{p-p{to)-,p-p{to)-) = 9{p{to)-,p{to)-) + g{p,p-2p{to)-) 

= 9ip{to)-,p{to)-) + g (^,p-2 ^hm p(to) 

Now, note that the covector p = ^limj_^^+ 'Pq(t)^ p(to)- ^ ^q(to) ^'^ch that the expression 
g{p — p{to) ^ , p — p{tQ) ^) is minimized among all the covectors belonging to D~^{t^^y Therefore, 
the derived jump rule (|T3p can be stated as follows: the "post-impact" momenta p(to)+ 
is such that the kinetic energy corresponding to the difference of the "pre-impact" and 
"post-impact" momentum is minimized among all the covectors satisfying the constraints. 
This is an appropriate version for generalized constraints of the well-known jump rule for 
perfectly inelastic collisions . This is even more clear in the holonomic case, as is shown 
in Subsection 4.2. 
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Remark 4.4 So far, we have been dealing with impulsive constraints. More generally, 
we can consider the presence of external impulsive forces associated to external inputs or 
controls. Then, equation (0) must be modified as follows 



{PA{to)+ - PA{to)- - P'a{U)) dql^t^) e lim = 



(16) 



where P^(to), ^ A < n, are the external impulses at time to. Observe that if g(to) is a 
regular point then 

PA{to)+ = PA{to)- + PA(io) - C''^jBg''''u;,APcito) , A=l,...,n, (17) 

i,j,A,B 



and, if q{to) is a singular point, we have 



PA{to 



p^(to)-+PA(io)-lim ( ^ 



lit) 



iPcito)-+Pcito)), A = l,...,n. 

(18) 



4.2 The holonomic case 

We show in this section a meaningful interpretation of the proposed jump rule (0) in case 
the codistribution D is partially integrable. 

Let us consider a trajectory q{t) G Q which reaches a singular point g(to) ^ "S* and falls in 
either Case 1 or Case 3. Since Q = R, we have that q'(to) ^ L, where L is the leaf of D which 
contains the regular posterior points of the trajectory q{t). On leaving q{to), we have seen 
that the trajectory suffers a finite jump in its momentum in order to satisfy the constraints 
imposed by D, which in this case implies that the trajectory after time to belongs to the leaf 
L. Consequently, the jump can be interpreted as a perfectly inelastic collision against the 
"wall" represented by the leaf L\ 



Let us see it revisiting Example 4.1. 



Example 4.5 Consider again the situation in Example [4.1| . If the motion of the particle 
starts on the left half-plane going towards the right one, then it is easy to see that = {0} 

and D^y) = spa.ia{dx — dy} . As D^oy) ^ -^(oy)' a jump of momenta is possible. In fact, if the 
"pre-impact" velocity (ioj^o) does not satisfy xq = yo, the jump occurs and is determined 
by Av{to) G D^oy) ^(^o) ~ vi^)- Consequently, we obtain 



y{4 



2 
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We would have obtained the same result if we had considered that our particle hits, in a per- 
fectly inelastic collision, against the "wall" represented by the half-line of slope 1 contained 
in {x > 0} passing through the point {0,y). 

If the particle starts on the right half-plane towards the left one, the roles are reversed and 

^{Oy) ~ ^' "^(Oj/) ~ 'f^J'' have that D'^^y-^ ^ ^{Oy) therefore there is no 

jump. 



5 Examples 



Next, we are going to develop two examples illustrating the above discussion. First, we treat 



a variation of the classical example of the rolling sphere pTl , Secondly, we take one 
example from Chen, Wang, Chu and Chou M. 



5.1 The rolling sphere 

Consider a homogeneous sphere rolling on a plane. The configuration space is Q = x 
5*0(3): (x, y) denotes the position of the center of the sphere and {ip, 9, ip) denote the Eulerian 
angles. 

Let us suppose that the plane is smooth if a; < and absolutely rough if x > 0. On the 
smooth part, we assume that the motion of the ball is free, that is, the sphere can slip. But 
if it reaches the rough half-plane, the sphere begins rolling without slipping, because of the 
presence of the constraints imposed by the roughness. We are interested in knowing the 
trajectories of the sphere and, in particular, the possible changes in its dynamics because of 
the crossing from one half-plane to the other. 

The kinetic energy of the sphere is 

where uj^, ^y and are the angular velocities with respect to the inertial frame, given by 

Ux = cos ip + if sin 6 sin , 
Uy = Osmip — ip sin cos ip , 
ujz = ^ cos 6 + ifj . 

The potential energy is not considered here since it is constant. 

The condition of rolling without sliding of the sphere when x > implies that the point of 
contact of the sphere and the plane has zero velocity 

(j)^ = X — rujy = , 
0^ = y + ru;x = 0, 
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where r is the radius of the sphere. 



Following the classical procedure, we introduce quasi-coordinates "gi", "^2" and "gs" such 
that "gi"= ujx^ "g2"= '^y and "^3"= oJz. These last expressions only have a symbolic meaning 



where we interpret "(igj" and "7^", 1 < i < 3, as adequate combinations of the differentials 



and partial derivatives, respectively, of the eulerian angles. Note that {^,^5^, 

"^"} and {dx, dy, dz, '^dq^" , '^dq'^" , ''dq^^^ } are dual bases. Moreover, observe that, 
from the nonintegrability of the constraints, the differential forms "(ig"' do not represent 
exact differentials. 

The non-holonomic generalized differentiable codistribution D is given by 



D(x,y,(f,,e,tp) 



{0} , if X < , 

span {dx — rdq'^jdy + rdq^} , if x>0. 



The intersection of the regular set of the generalized codistribution and the (x, ?/)-plane has 
two connected components, the half-planes Ri = {x < 0} and R2 = {x > 0}. The line 
{x = 0} belongs to the singular set of D. 

On Ri the codistribution is zero, so the motion equations are 

mx = , 
my = , 



mk^uj^ = 0, (19) 
mk'^Ljy = , 
mk'^u)^ = . 



On i?2 we have to take into account the constraints to obtain the following equations of 
motion 



mx 


= Ai, 


my 


= A2, 


mk'^Ux 


= r\2 , 


mk'^tOy 


= — rA 


mk'^LUz 


= 0, 



(20) 



with the constraint equations x — rujy = and y + rujx = 0. One can compute the Lagrange 
multipliers by an algebraic procedure described in ]TB|. 



Suppose that the sphere starts its motion at a point of Ri with the following initial conditions 
at time t = 0: xq < 0, yo, xq > 0, yo, {uJx)o, (t^y)o and {uJz)o- Integrating equations (pJ]) we 
have that if x{t) < 

x(t) = Xot + Xo, 
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UJy{t) 



mt + yo , 

(^'.t)o • 



(21) 



At time t = —xq/xq the sphere finds the rough surface of the plane, where the codistribution 
is no longer zero and it is suddenly forced to roll without sliding (Case 1). Following the 
discussion in Section 4, we calculate the instantaneous change of velocity (momentum) at 
a; = 0. 

First of all we compute the matrix C 



1 -r 
1 r 



/ 1 











\ 




1 


0\ 





1 
















1 








k-' 













r 


















— r 
















k-^ J 


\ 





J 



{1 + r'k-') 



1 
1 



Next, a direct computation shows that the projector V does not depend on the base point 





r2 














rk'^ 




r2+fe2 








r2+fc2 
-rfc2 








' r'^+k'^ 






Therefore, we have 



(Px)- 
(Py)- 
iPi)- 

{P2)- 

iPs)- 



r 


o\ 














r2+fc2 








1/ 








_|_ ^2 ' 

^^(Py)o ~ ^(Pi)o 

^2 _|_ ^2 ' 

-r/c^(py)o + /;;^(pi)o 

^2 _|_ ^2 ' 

rF(;;,)o + A--(y^2)o 



y.2 _|_ ^2 



= (P3)o ■ 



Now, using the relation between the momenta and the quasi-velocities 

Px ^i:Py = = k^U;x:P2 = k'^UJy.Ps = k^i^z : 
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we deduce that 



y+ 



^2 _j_ ' 
_j_ ^2 ' 



= ^^^^^I^^. (22) 



^2 _|_ ^2 ' 

Finally, integrating equations (pOD at time t = —xq/xq with initial conditions given by (|2 
we obtain that if t > t 

_l_ 



rip + P(u;^)o 

J.2 _|_ ^2 



5.2 Particle with constraint 

Let us consider the motion of a particle of mass 1 in subjected to the following constraint 

= (y^ — — z)x + {z — y'^ — xy)y + xi = . 

In addition, let us assume that there is a central force system centered at the point (0, 0, 1) 
with force field given by 

F = —xdx — ydy + (1 — z)dz . 
Then, the Lagrangian function of the particle is 

L = T -V = ^{x^ + f + z"" + x' + y^ + z^ -2z) , 

and the constraint defines a generalized differentiable codistribution D, whose singular set 

is S = {{x,y,z) : X = 0,z = y'^}. 

On R, the regular set of D, the dynamics can be computed following the standard symplectic 
procedure to obtain Tld = Tl + \Z, where 

^ d d d d d . . d 

= +?/^ + -3;— -y— - (2;- 1)— , 

ox oy oz ox oy oz 
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and A is given by 

^ Tl^cJ)) —2xx^ + yijx — 2yy'^ — xy'^ + yz + x^ + y^ — yz + x 
Z{(f)) — x'^ — zy + {z — y'^ — xyY + x"^ 

Consequently, the motion equations on R are 

mx + X 
my + y 
mz + z — 1 

with the constraint equation = 0. 

From the discussions of 0, we know that in this case there is an integral surface, C, of the 
constraint 0, that is, a surface on which all motions satisfy the constraint. This surface is 

C = {{x,y,x) : z — x^ — y"^ + xy = . 

Note that S G C. Therefore, if a motion takes place on the cone-like surface C, it is confined 
to stay on this critical surface, unless it reaches a singular point. In this case, the space of 
allowable motions is suddenly increased (in fact, TM^), and the motion can "escape" from 
C. In addition, this proves that the unique way to pass from one point of the exterior of the 
C to the interior, or viceversa, is through the singular set S. 

In particular, we are interested in knowing 

(i) Is there any trajectory satisfying equations ( pSf ) which passes through the singular set? 

(ii) if so, which are the possible momentum jumps due to the changes in the rank of the 
codistribution D7 



= - X - 2;) , 

= Xiz~y^-xy), (25) 
= Aa; , 



So far, we do not know an answer for the question of the existence of a motion of ([25| ) 
crossing S. It seems that on approaching a singular point, the constraint force can become 
increasingly higher (p^). Consequently, this force possibly "disarranges" the approaching of 
the motion to 5*. Numerical simulations are quite useless in this task, because of the special 
nature of the problem: the hard restriction given by the fact that a motion crossing the 
cone-like surface C must do it through the singular part S. Indeed, the numerical simulation 
performed in crosses the surface C through points which are not in S. 

Concerning the second question, let us suppose that there is a trajectory of the dynamical 
system (^), q{t) = {x{t),y{t), z(t)), that passes through a singular point at time to, i.e. 
x(to) = and z(to) = y'^ito). The rank of the codistribution D at the immediately preceding 
and posterior points is 1, meanwhile at g(to) it is (Case 3). So, a possible jump of the 
momentum can be induced by the change in the rank of D. 

A direct computation shows that the projector V depends exphcitly on the base point q E Q. 
Equivalently, we have that D^^it^) depends strongly on the trajectory q{t). In fact, taking 
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two curves qi{t), q2{t) passing through q'(to) at time t^, and satisfying Xi{t) ^ Zi{t) — yf{t) 
and Z2{t) — ^ X2{t) when t ^ respectively, one can easily see that -D^^^^^ 7^ ^q2{to) 
(the expression /(t) <^ fl'(^) when t — > means that limj^j+ f{t)lg{t) — 0). 

Consequently, we are not able to give an answer to question (ii) (in case the first one was 
true) unless we assume some additional information: for example, that the balance between 
x{t) and z{t) — y'^it) is the same for t and t — > t^. In such a case, -D"^^^^ = -D^^^^^-j and we 
would conclude that there is no jump. In mechanical phenomenae of the type sliding-roUing, 
as the ones studied in Section 5.1, this kind of "indeterminacy" will not occur in general. 

In spite of the fact that the most natural thing in this case seems to be to think that there 
is no jump of momenta, a mathematical explanation of it is still to be found. 
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